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ABSTRACT 


In experimental-numerical technique of the present study, the aim is to find the initia- 
tion and propagation toughness of fast moving interlaminar cracks in GFRP laminates. 
Experimental technique consists of a specimen, the front face of which is bonded to a 
rigid block and the cantilever plate on the rear is screwed to a load bar. A striker bar 
impacts the load bar, which in turn impacts the specimen and causes the crack prop- 
agation. Strain gauges mounted on the load bar monitor the end displacement of the 
load bar. Since the specimen is screwed to the load bar, the end displacement of the 
load bar is same as the deflection of the cantilever end. Two strain gauges bonded ahead 
of the crack tip on the specimen monitors the crack velocity. Thus the experimentation 
provides the deflection of the cantilever end, the crack velocity and the crack initiation 
time. 

This data along with the material properties are used as input to the FE code to sim- 
ulate the crack propagation. Gradual nodal release method is used to model the crack 
propagation. The variation of J integral with time is obtained from the program. It’s 
value steadily increases till the initiation time. Value of J integral at the initiation time 
is treated as the initiation toughness (Ji„i). At the initiation time the crack propagation 
algorithm is called. The J integral suddenly drops at that point and stabilizes after some 
time. This stabilized of J is treated as propagation toughness Jprop- The initiation and 
propagation toughness of a crack moving at a speed of 550-950 m/s were found to be 
140-350 J / rr? and 30-95 J / vn? respectively. These values are very much lower than the 
quasistatic interlaminar toughness ( 1050 J/m^). 
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Chapter 1 
Introduction 


1.1 Introduction 

The word "composite" in composite material signifies that two or more materials are 
combined on a macroscopic scale to form a useful material. The key is the macroscopic 
examination of a material. Different materials can be combined on a noicroscopic scale, 
such as in alloying, but the resulting material is macroscopically homogeneous. The 
advantage of composites is that they usually exhibit the best qualities of their constituents 
and often some qualities that neither constituent possesses. The properties that can be 
improved by forming a composite material include 

• Strength 

• Stiffness 

• Weight 

• Fatigue life, etc. 

Therefore the fibre reinforced composite materials have seen extensive application espe- 
cially in aerospace industry, where the weight saving is most essential. They have good 
strength, stiffness and fatigue life but are suceptible to delamination when a foreign ob- 
ject impacts it at high velocity. On a metallic sheet (e.g. steel, aluminium) only a dent 
is formed when it is impacted by a foreign object with low or intermediate velocity, the 
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damage is local. Where as in case of composite materials it causes delamination. Sev- 
eral such impacts leads to extensive delamination and then the material integrity is in 
question. Therefore it is of our interest to study the dynamic fracture behaviour of the 
composite under impaet loading. 

Dynamic fracture phenomenon has several important featmes. The mathematical models 
are complex than the static models. The boundaries of the body changes with time. From 
experimental point of view, many parameters have to be measured accurately with in a 
very short period. 

1.2 Literature Survey 

A large number of investigations have been carried out on dynamic crack propagation 
employing through the thickness crack in a large plate. Stress waves are continuously 
generated at the crack tip of a fast moving crack. As these stress waves reach the bound- 
ary of the specimen, they are reflected. Some of these reflected stress wave come back 
to the region just in front of crack tip. Experiments with large size plate are controlled 
such that the monitoring of important parameters is complete before the reflected waves 
arrive. In other words, the experiments are controlled to isolate the effect of specimen 
edges. In case of an interlaminar crack in a FRP laminate, which is usually not thicker 
than 6 mm, the free surfaces axe always close to the crack tip and there is no way one 
can isolate the region of interest from the reflected waves. Therefore a new approach is 
required to find dynamic toughness of interlaminar cracks.In this section, investigations 
on the large plates with through the thickness axe briefly reviewed. 

Initial and extensive investigations on crack propagation of through the thickness crack in 
a large plate were done by several groups. Ravi Chandar and Knauss (1982, 1984a, 1984b, 
1984c, 1984d) have done an exhaustive study of dynamic crack propagation phenomena 
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using the method of caustics for Homalite-100. This included study of crack initiation and 
arrest, microstructural aspects, crack branching and interaction of stress waves with the 
crack tip. Rosakis, Duffy and Pruend (1984) performed dynamic crack propagation ex- 
periments on double cantilever beam specimen using wedge loading. Zehnder and Rosakis 
(1990) used optical method of reflected caustics combined with high speed photography 
to investigate the dynamic fracture initiation and propagation in 4340 steel specimen. 
These investigations, based on caustic have been very effectively used to study the frac- 
ture phenomena in a large plate with through the thickness crack. The method is not 
likely to be effective for the interlaminar crack in a slender sheet because the strength of 
the interlaminar bond is very weak and it will introduce very little strain in the bonded 
sheet making the size of the caustic negligibly small. 

Ravichandran and Clifton (1989) developed a special technique to study the initiation 
and propagation of crack in the steel under dynamic impact loading. They presented a 
plate impact experiment and an associated finite difference model to study the fracture 
process that occurs in sub-micron loading. A disc containing a prefatigued edge crack 
on the midplane up to half way across the diameter is impacted by a thin flyer plate of 
same material. A compressive pulse propagates through the specimen and reflects from 
the rear siurface as a tensile pulse of 1/xs duration. The motion of the rear surface is 
monitored by using the laser interferometer. 

Sim and Grandy (1988) investigated dynamic delamination fracture toughness in a [90/0]5s 
T900/934 graphite/epoxy laminate using impact loading. Delamination cracks of differ- 
ent sizes were embedded at the midplane of the composite specimen. The threshold 
impact velocity that causes propagation of delamination crack was used in the dynamic 
analysis with the finite element method. From the finite element solution, the time- 
history of the strain energy release rate was calculated. The critical strain energy release 
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rate was taken equal to that of maximum value of the response history. 

Berger and Dally (1990a) used a series of strain gauges ahead of the crack tip to monitor 
the strain and crack propagation. Berger, Dally and Sanford (1990b) also used strain 
gauges ahead of the crack tip to determine dynamic stress intensity factor associated 
with a propagating crack. 

Takeda et.al., (1982) used high speed photography to measure the speed of interlaminar 
in composite laminates. They observed that the interlaminar crack in composite lami- 
nates moves at 200-500 m/s, when impacted by foreign bodies. Freund (1990) found that 
the crack velocity is subsonic i.e. less than Rayleigh wave speed. 

Nishioka and Atluri (1983) studied the use of path independent J-integral for dynamic 
crack propagation by the finite element method. Other path integrals were also investi- 
gated along with J-integral. Numerical results showed that combined use of J-integral 
and the finite element method is a useful tool to obtain the firacture parameters such as 
stress intensity factors and energy release rates. 

Kolednik (1991) presented his theoretical study for physical interpretation of the J-a 
curves for elastic-plastic fractmre. He derived the difference between the exact initiation 
toughness and crack growth toughness for a very low velocity crack using energy'" balance 
under quasi-static conditions. He considered three point bend specimens each of which 
consists of two parts glued together along the ligament. The analysis was made on large 
specimens. 

Verma et al. (1995) investigated SIF in a DCB specimen made of thin cantilevers of steel 
bonded with epoxy through a combined scheme of measuring the strain near the crack 
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tip and analyzing the experimental data using FEM. In the finite element programme 
he developed a relationship between stress intensity factor and strain near the crack tip. 
The strain measured near crack tip with a strain gauge is used in the FEM programme 
to obtain the SIF of the DCB specimen. Lovi (1993) improved this method by using 
two strain gauges on each cantilever to neutralize the possibility of recording the bending 
strains in the cantilevers which might develop during the specimen preparation. 

Verma (1995) developed a combined experimental and numerical technique to determine 
interlaminar dynamic firacture toughness under impact loading. An interlaminar crack 
is propagated at very high speed in DCB specimen made of two steel plates which are 
bonded together by epoxy with a precrack. Strain gauges axe bonded ahead of crack tip 
to monitor the crack velocity. The experimentally obtained data such as the end deflec- 
tion of the cantilever, crack velocity and crack initiation time are used as input to the 
FE code to simulate the crack propagation, and to find the initiation and propagation 
toughness. 

Ramakrishna (1997) applied the technique developed by Verma to composite laminates. 
He obtained initiation and propagation toughness of GFRP laminates which are of rela- 
tively poor quality. Babu(1998) further exploited the same technique to find the initiation 
and propagation toughness of GFRP laminates of relatively superior quality. 

Raman P. singh et al. (1996) described various experimental observations for dynamic 
intersonic decohesion of bimaterial interfaces. Two separate but complementary optical 
methods are used in conjunction with high speed photography to explore the nature of 
the large scale contact and shock wave formation at the vicinity of running cracks in two 
different bimaterial systems. 
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R.W.Truss et al. (1997) had attempted to find interlaminar and intralaminar fracture 
toughness of uniaxial continuous and discontinuous carbon fiber/epoxy composites, using 
compact tension and double cantilever beam test geometries. The discontinuous carbon 
fiber /epoxy composites have been found to have a slightly misalignment of the fibers 
from the average fiber direction and this misalignment was found to increase both initi- 
ation fracture toughness and to greater extent the propagation fracture toughness. The 
increase in fracture toughness in discontinuous carbon fibers/epoxy samples was due to 
fibers bridging the crack and this has been modelled as if the the bridging fibers provide 
an increase in compressive stress across the crack. 

A.Pegoretti et al. (1998) investigated effects of the temperature and fiber surface finish- 
ing on dynamic interfacial fracture toughness of nylon-6/E glass fibers composites. They 
have conducted experiments on embedded single fiber micro composites. For the unsized 
E-glass fibers an interfacial toughness of 165 J/m^, while for the polyamide and epoxy 
compatible E-glass fibers debonding energy value was found as 340 J/m^. Both types of 
E-glass fiber micro composites showed a strong decrease of interfacial toughness as the 
temperature increased. 

John Lambros and Rosakis A.J (1997a) investigated dynamic delamination of thick fiber 
reinforced polymeric matrix composite laminates using optical techniques and high speed 
photography. They had used 65% fiber volume fraction of Graphite/epoxy laminates 
consisting of 48 plies. Square plates of 152mT7r x 152mm dimensions were impacted in 
an out-of-plane configuration using high speed gas gun. Real time imaging of the lami- 
nate out-of-plane displacement was performed using lateral shearing of interferometer of 
Coherent Gradient Sensing in conjunction with high speed photography. Delamination 
speeds up to 1800 m./s were observed. 
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John Lambros and Rosahis A. J (1997b) investigated dynamic crack initiation and growth 
in unidirectional graphite/epoxy plates with high speed photography. Edge notch plates 
are impact loaded in one point bend configuration using drop weight tower. Initia- 
tion fracture toughness data were reported and significant dynamic effects are observed 
through emission of stress waves from the propagating crack tip. 

1.3 Present Work 

Verma(1995) developed a technique to determine the interlaminar toughness of two steel 
specimens bonded with epoxy. The technique employs a combined experimental and nu- 
merical work. Ramakrishna (1997) , Babu(1998) developed the same technique to GFRP 
laminates. In this work the technique has been further exploited to determine the dy- 
namic initiation toughness and propagation toughness. 

The technique consists of applying an impact load on a specimen with a precrack (Fig 
1.1). The portion that is on the front side of the precrack remains straight during exper- 
imentation because it is bonded to a solid block. The rear rear portion of the specimen 
works as a cantilever plate. When an impact load is applied on the cantilever plate, 
the crack propagates. Strain gauges are mounted ahead of the crack tip to monitor the 
response. The experimental work determines the end deflection of the catilevr, crack 
velocity and the crack initiation time. These are used as input to the FE code developed 
by Verma (1995) to simulate the crack propagation. 

Chapter 2 of the thesis describes the experimental setup and the details of the specimen. 
Chapter 3 describes the brief outline of the FE code developed by Verma (1995). Results 
are discussed in chapter 4 and finally the conclusions and the scope for further work are 
stated in chapter 5. 
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Chapter 2 

Experimental Technique 


2.1 Introduction 

In the experimental - numerical technique of the present study, the experiments are con- 
ducted to determine the deflection of the cantilever end, crack propagation velocity and 
initiation time under impact load on the specimen. (Fig 1.1). This chapter deals with 
the specimen details and its preparation, experimental setup, experimental measurements 
and data analysis. 

The technique developed by Verma(1995) is adopted to the present work. In this tech- 
nique a stress pulse is generated in a cylindrical load bar (Fig. 2.1) by impacting it 
with a striker bar, accelerated in an air-gim. This load bar impacts the cantilever end of 
the specimen, whose other end is bonded to a very high inertia solid block. Two strain 
gauges are mounted in symmetrically opposite direction, on the surface of the load bar to 
monitor the stress pulses that is acting on the specimen. Two strain gauges are mounted 
ahead of the crack tip, to measure the crack velocity. Oscilloscope is used to record 
the data from the strain gauges. All the strain gauges are connected to the oscilloscope 
through the bridge circuits. The details of the experimental setup would be presented in 
Sec 2.3. 
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2.2 Specimen Details 


The specimen used in the present work is unidirectional glass fibre reinforced epoxy com- 
posite. This section describes the basic raw materials used for the specimen, preparation 
of the laminate, geometry of the specimen, and preparation of the specimen. 

2.2.1 Preparation of the Laminate 

The basic raw materials for the preparation of the laminate are glass fibre and epoxy 
resin mixture. The composition of the epoxy resin is as follows: 


Araldite 

LY556 

100 parts by weight 

Hardener 

HT976 

35 parts by weight 

Accelerator 

XY73 

1 part by weight 


Coupling agent ( 7 - Ami no propyl triethoxy silane ) 0.5 parts by weight 

Unidirectional glass fibres ( of diameter 5-20 /im ) reinforce the epoxy resin matrix. 
The glass fibre unidirectional prepegs are obtained in the semicured state firom the uni- 
directional prepeg machine available in the laboratory. Number of lamina can be stacked 
together to get the laminate of required size. In the present work 52 of such lamina (of 
length 200 mm) are stacked with a precrack ( thin tefion insert, 0.02 mm thick) placed 
exactly at the mid plane and to one end of the laminate. This laminate is then cured 
in a hydraulic press. The laminate is placed between the two platens of the press and is 
heated to 120 ° C. Then the pressure is gradually raised to 0.7 Mpa (~ 7 ata) in half an 
hour. At the same pressure and temperature it is allowed to cure for four hours. Then 
keeping the pressure same, the temperature is increased to 150*^ C and is maintained for 
3 hours. It is then allowed to cool to room temperature. The cured laminate is taken out 
of the press and is cut to required sized specimens by a diamond cutter. 
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2.2.2 Specimen Geometry 

The geometry of the specimen is as shown in Fig.2.2. Figure 2.3 shows the isometric 
view of the same. It is 25 mm wide and 180-200 mm long with an edge precrack on the 
mi d plane of the specimen. The portion that is on the front side of the precrack remains 
straight during experimentation because it is bonded to the solid block. The rear portion 
of the specimen works as a cantilever plate. The thickness of the cantilever is about 4 
mm. A hole of (f) 6mm is drilled on the cantilever end in order to fix it to the load bar. 
The center of this hole is the loading point on the cantilever. 

2.2.3 Crack Sharpening 

The as-cast precrack in the specimen is not very sharp due to the finite thickness of the 
teflon insert. So further extension of the crack is needed to sharpen the tip. It is done 
with the help of a fixture and pressure plates as shown in Fig 2.4. Two pressure plates are 
cut to 2 mm depth leaving a width of 8 mm at each end, to ensmre that high compressive 
stress exists at a few millimeters away from the as-cast tip. Before the presstue plates 
are tightened between the two jaws of a vice, they are placed on their side face over the 
third base plate whose top surface is also grounded fiat. This ensures that the edges of 
the pressure plates are parallel to the crack front. The specimen is placed between the 
pressure plates in such a way that the tip of the pre crack is out by 2-3 mm beyond 
the edges of the pressure plates. The specimen between the pressure plates is pressed 
properly between the jaws of the vice. Now with a cautious pressing of a sharp chisel at 
the mouth of the precrack the crack is extended up to the edges of the pressure plates. 
The tip of the extended crack is sharp. 

2.3 Experimental Setup 

Figure 2.5 shows the photograph of the overall experimental setup, the schematic diagram 
of which is as shown in Fig. 2.1. The firont plate of the specimen is bonded to a rigid block, 
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which in turn is clamped to the rigid base using C-clamps. The cantilever plate of the 
specimen is screwed to the load bar. Two strain gauges are bonded at diametrically op- 
posite locations of the load bar. These are connected to an oscilloscope through a bridge 
circuit. Each one of the velocity strain gauges which are located ahead of the crack tip 
on the specimen is also connected to the same oscilloscope through another bridge circuit. 

The striker bar accelerated in the barrel of the air-gun impacts the load bar. This pro- 
duces a compressive stress pulse in the load bar which propagates towards the specimen. 
When the load bar impacts the specimen some part of its energy is transmitted to the 
specimen and the rest is reflected back as tensile pulse. The energy that is transmitted 
to the specimen causes the deflection of the cantilever plate and in turn propagates the 
crack. The strain gauges mounted on the load bar records the incident and reflected 
pulses through the oscilloscope. The velocity strain gauges on the specimen, which are 
connected to the oscilloscope, record the crack propagation history. 

2.3.1 Solid Block 

As shown in Fig.2.1 the front plate of the specimen is bonded to a rigid block of Mild 
steel of dimensions Ybrnm x 75Tnm x ISOmm. The block is placed on a rigid base plate; 
after aligning the specimen to the load bar, the block is firmly clamped to the base plate 
with the help of C-clamps. Care is to be taken to make sure that the crack plane is 
exactly perpendicular to the load bar; otherwise, the specunen will be preloaded and the 
numerical simulation would be inaccurate. 

2.3.2 Striker and the Load bco* 

Applying an impact load on the specunen includes an air-gun, a striker bar and a load bar. 
The striker and load bar are made of cold rolled mild steel of diameter 19 mm. The load 
bar is aligned properly with the center line of the barrel of the air-gun. The striker bar is 
accelerated in the air-gun and its impacting face is made spherical (radius ~90 mm) . It 
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impacts the load bar, which in turn impacts the specimen. When the striker bar hits the 
load bar, a compressive stress pulse is produced in the load bar which travels towards the 
specimen. A part of this compressive stress pulse is transmitted to the specimen and the 
rest is reflected back as a tensile pulse. Two strain gauges (120f2, 6 mm gauge ) mounted 
at diametrically opposite locations on the load bar are used to monitor the stress pulses. 
These strain gauges, connected to the oscilloscope through a bridge circuit monitors the 
incident and reflected stress pulses. 

2.3.3 Bridge Circuit 

The purpose of bridge circuit is to convert the resistance change in strain gauges into a 
voltage difference. This voltage difference can be recorded in the oscilloscope. Figure 2.6 
shows the details of the bridge circmt used. There are four arms in each bridge. One 
strain gauges is connected to each of these arms. They are named R], R 2 , R 3 , R 4 . In an 
half bridge circuit two strain gauges are active gauges, (i.e Rj and R 3 ) and in a quarter 
bridge circuit only one gauge (i.e Ri) is active. The remaining gauges are dummy gauges. 
These dummy gauges are provided by mounting them on a plate of the same material as 
that of the load bar, which acts as a heat sink for the heat generated in the gauges when 
current passes through the gauges. The half bridge circuit is used to monitor the stress 
pulses in the load bar. In this case the two active gauges are the two gauges mounted on 
the load bar at diametrically opposite locations. 

The quarter bridge circuits are used to monitor the crack propagation history through 
velocity strain gauges. Since there are two velocity strain gauges bonded on the specimen 
two quarter bridges are needed to monitor the crack propagation history. 

The bridge circuit is balanced to zero volatage output by connecting 1 resistance for 
coarse adjustment and 1 M O resistance for fine adjustment, as shown in the circuit. To 
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calibrate the circuit a calibration resistance ( R c = 47.0 KO) is connected in parallel to 
any active gauge through a switch. The equations involving the calibration of the bridge 
circuit are as follows. 

For a balanced bridge circuit the output voltage Ae is given by 


Ae 


R^.R2 ARi AR 2 ARs _ AR4 
(R1+R2)' -^2 Rz Ra 


( 2 . 1 ) 


where E is the input voltage. It can be seen from the equation that similar ( both positive 
or both negative ) changes in resistance of opposite arms of the bridge circuit are added 
and dissimilar ( one positive and other negative ) changes are cancelled. Thus by taking 
the active strain gauges at the opposite arms, only compressive pulse is recorded in the 
load bar and bending effect is nullified, if there is any. 


The relation between the strain in the strain gauge and corresponding change in its re- 
sistance is governed by the following equation 


c _ 

^ AL/L 


( 2 . 2 ) 


where Sg is the gauge factor and AR/ R is the strain recorded by the strain gauge. 
Rearranging the above equation. 


AL _ AR/R 

L Sg 
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For Half Bridge 


The change in resistance of the arm AB after connecting Rc is given by 

This implies 

m 

Corresponding to this change in resistance a voltage difference ( Calibration 
Voltage, Vc) develops between terminals A and C. Similar resistance change will also 
occur in Ri and R3, when the load bar experiences the load pulses. From Eqs. 2.3 and 
2.5 the calibration voltage corresponds to the strain 


Cc = AI/T 


2(i?i + Rc)Sg 


( 2 . 6 ) 


The strain value, ec will thus correspond to the calibration voltage, Vj.. A factor of 1/2 
is introduced to take the average of the strains recorded by the strain gauges Ri and R3 
which are at the opposite arms of the bridge circuit. By the linear relationship between 
the voltage drop across AC and strain in the strain gauges bonded on the load bar, the 
strain in the load bar corresponding to a voltage V recorded on the oscilloscope can be 
given as 


€ = 



Ri V 

2 (i?i + Rc) Sg Vc 


(2.7) 


For Quarter Bridge 

The analysis for this circuit is also same as that of the the half-bridge circuit except that 
there is only one active strain gauge (i?i), Rz acting as a dummy gauge in this case. Since 
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there is only one active gauge €c is 


AL/L 




{R,+R,)S, 


and strain at the crack velocity strain gauge is 


. = = 


Ri V 


{R, + R,) Sg K ■ 


( 2 . 8 ) 


(2.9) 


The strain gauges are connected to the Wheatstone bridge through a specially designed 
cap and stud arrangements to have press contact at the terminals ( Verma (1995) ). 
Coaxial wires are used to connect the circuit to all the gauges to shield the fast changii g 
signals. A bridge is enclosed within a conducting box of aluminium to have shielding from 
spurious signals. The groimds of all coaxial cables and the oscilloscope are connected to 
the aluminium enclosures to have the same common ground of the entire measuring 
system. 

2.3.4 Oscilloscope 

Each one of the strain gauges is connected to the oscilloscope through the bridge cir- 
cuits. The oscilloscope records the resistance change in the strain gauges as the potential 
difference. The oscilloscope used in the present work is a 4-charmel Digital Storage Oscil- 
loscope ( Model 1624, Gould Inc., U.K). It has 8-bit resolution in vertical direction and 
16-bit resolution in horizontal direction. The sensitivity of the Oscilloscope is 0.025 mV 
on vertical axis and 0.25 /xs on horizontal axis. It is integrated with a personal computer. 
A PC compatible software is used to transfer the data from the oscilloscope to the PC 
for further analysis. 

2.3.5 Deflection of Cantilever end 

Figure 2.7 shows a typical record of the incident and reflected pulses in the load bar. From 
this experimental record the deflection of the cantilever end can be obtained using one 
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dimensional wave propagation theory. Figure 2.8 shows the time-distance (t-x) diagram 
for the propagation of stress pulses in the load bar. When the striker bar impacts the load 
bar, the compressive incident stress pulse (crj) propagates towards the specimen and is 
recorded at the location 1. The reflected tensile stress pulse (<73) is recorded at the same 
location but at a diferent time and is depicted by point 3 in the t-x diagram. The aim is 
to express particle velocity V2, at point 2 in terms of experimentally recorded cri and 0-3. 
One dimensional wave equation along the characteristic with the acoustic impedance p c 
is expressed as 

dcr-pcdv=0 ( along -1-ve characterstic) 

dcr-l-pcdv=0 ( along -ve characterstic) 

where a is stress, v the particle velocity, c the longitudinal wave velocity, and p the 
density of the load bar material. 

Using the above concept the relation along characteristic 1-2 can be obtained as 

<72 — p C t)2 = <7i — p c Vi . (2-10) 

Along the characteristic 1-5 the relation is 

ai+pcvi = ar^ + pcvz . ( 2 . 11 ) 

But <75=0 and ^5=0, because the load bar is initially at rest and stress wave never reaches 
point 5. Then the above equation becomes 

<7] = — p c Vi . (2-12) 

Substituting Eq.2.12 in Eq.2.10, 

<72 — p c «2 = 2ai . (2.13) 
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Relations along the characteristics 2-3 and 3-4 are 


a 2 + p c V 2 = (73 + p c Vs 


(2.14) 


and 


(7z - p C Vs = (74 - p C Vi . 


(2.15) 


It is worth noting that <74 = 0 and 1)4 = 0 at point 4, because the striker bar and the load 
bar are made of the same material and diameter and the striker comes to rest. Then the 
above two equations yield 

a2 + pcv2 = 2(73 (2.16) 


Adding equation 2.13 and Eq.2.16, one obtains 


(72 = CTj -h (73 


(2.17) 


By substituting <73 from the above equation in Eq. 2.13, 


<73 — <7i 
= 


( 2 . 18 ) 


fi c 

Since is compressive and crs is tensile in nature the particle velocity of cantilever end 
is found by taking the sum of the absolute value of the incident and the reflected pulses. 
The displacement U 2 of the load bar end can be found by integrating the velocity. It is 
given by 


U2{t) = f V2 dt = f {- —) (2.19) 

Jo Jo P c 

Since the specimen is screwed to the load bar the displacement of the load bar end 
becomes the defelction of the cantilever end. Figure 2,9 shows a typical variation of the 
deflection of the cantilever end with time. 
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2.3.6 Crack Velocity 

Verma(1995) developed a technique to find the dynamic interlaminar toughness of steel 
DCB specimen bonded with epoxy; luider impact loading. He used strain gauges bonded 
ahead of the crack tip to monitor the crack propagation history. The strain gauges used 
in his work were supplied by Tokyo Sokki Kenkyujo Co. Ltd, Japan. They have a gauge 
length of 0.2 mm and a gauge factor of 2.05. The resistance of the gauges is 120 ±0.30. 
He placed the strain gauges at ±45°and at the mid plane of the cantilever. The technique 
worked quite reliably. 

Ramakrishna (1997) applied the same technique on GFRP laminates. It also worked reli- 
ably. However the interlaminar toughness of the GFRP of Ramakrishna’s work was poor 
due to the poor quality of the laminate. When Babu (1998) applied the same technique 
to GFRP laminates of substantially higher interlaminar toughness the technique based 
on more than 25 experiments, was found to be unreliable. Thus it was concluded that 
the technique developed by Verma, works only with specimen having low interlaminar 
toughness. 

Thus to overcome the problem a modified technique has been developed. In this tech- 
nique same strain gauges (0.2 mm gauge length, supplied by Tokyo Sokki Kenkyujo Co. 
Ltd, Japan) were used but they were bonded much closer to the crack plane and with 
a different orientation. The unnecessary portion of the strain gauge is cut and thrown 
away. ( Fig 2.10 shows the blown of views of the strain gauge prior to and after cutting 
). The strain gauge is then placed very close to the crack plane (Fig 2.11). Since the 
strain gauges are very close to the crack plane, the singular strain field gives a strain 
peak when the crack tip passes close to the strain gauge. There are two such velocity 
strain gauges the distance between the two strain gauges is measured accurately through 
travelling microscope; the velocity can be found ( Fig 2.12) using the formula 
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V elodty = — ^ 

h — ii 

where Cj is the location of strain gauge from the cantilever end and tj the time of the 
strain peak, monitored by the velocity stredn gauge. 

Verma (1995), using three velocity strain gauges, found that the crack velocity does not 
very much with time. Therefore in this study the crack velocity is assumed to be constant 
and only two strain gauges are employed. 

Bonding the strain gauges 

Initially the surface of the specimen is polished and degreased. Then the strain gauges 
are placed in their marked locations taking the help of a cello tape. Fewi-quick (Supphed 
by PediUte industries ) is used to bond the strain gauges. The strain gauges are placed 
very close to the crack plane and perpendicular to the crack plane. 

2.3.7 Initiation time 

The two strain gauges moxmted ahead of the crack tip feels the response of the crack 
thereby giving the crack velocity. The initiation time can be found by extrapolating 
these data. Figure 2.13 shows the details of the extrapolation. Extrapolation of the 
crack velocity to determine the initiation time can be justified as follows. 

John Lambros and Ares J. Rosakis (1997a) conducted an impact loading test on a one 
point bend (edge-notched unidirectional specimen with fibers parallel to the crack plane) 
configuration on unidirectional graphite/epoxy composite plates. For obtaining crack 
initiation and monitoring crack growth, an experimental technique with high speed pho- 
tography was used. The measured crack velocity, shown in Fig. 2.14, is in microsecond 
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domain; the crack, after initiation acquires high velocity immediately with rise time much 
less than 1 /xs. It is infered from the experimental result that the extrapolation invoked 
in this study in finding initiation time is not too bad. Further more it is clear from Fig 
2.14 that crack velocity does not change drastically after initiation time. It increases 
from 500-600 m/s in first 5 fis (20%). Thus if there is an error of 30% ( chosen higher to 
account for unkn own factors), in estimating the time taken by the tip of the precrack to 
the first strain gauge, the error in estimating the initiation time is only 1.10 /xs, for the 
crack to reach the first strain gauge at a distance of 3 mm with a crack velocity of 800 m/s. 


2.4 Closure 

The specimen details such as the basic raw materials, laminate preparation, crack sharp- 
ening and the overall experimental setup, which include rigid block, the stress pulses, the 
bridge circuit and the oscilloscope were discussed in this chapter. The experimentation 
provides 


• The deflection of the cantilever end 

• Crack velocity 

• Initiation time 


These are used as the input to the FE code ( developed by Verma (1995)). A brief outline 
of its formulation is presented in the next chapter. 


21 



SPECIMEN 



22 


Figure 2.1: Schematic diagram of the experimental Setup 









Fig 2.2 : Specimen Geometry 
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Figure 2.4: Crack sharpening fixture 
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Figure 2.5: Photograph of overall experimental setup 
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Figure 2.6: Details of the bridge circuit 
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Figure 2.7: Typical record of incident and reflected pulses in the load bar ( Expt.l) 
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Figure 2.8; Time - di.stanee ( t-x ) diagram 






(b) 

Figure 2.10: Blown up view of the strain gauges (a) before cutting , (b) after cutting 
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Figure 2.11: Position of the strain gauges with reference to the crack plane 
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Figure 2.12: Determination of crack velocity from the response times 
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Figure 2.13; Details of extrapolation to find initiation time corresponding to the length 
of the precrack (ao). 
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Chapter 3 

Numerical Analysis 


3.1 Introduction 

The data obtained from the experiments such as the initiation time, crack velocity and 
the cantilever end deflection are taken as input to the FE Code developed by Verma(1995) 
to simulate the crack under dynamic conditions. Verma developed the technique to find 
the interlaminar dynamic toughness of steel specimens under impact loading. Ramakr- 
ishna(1997) and Babu(1998) extended the same technique to the GFRP laminates. 

This chapter briefly describes the finite element formulation, the path independent inte- 
gral, the crack opening scheme and the flow chart of the code. 

3.2 Formulation 

The governing system of finite element equations for the linear dynamic response of an 
elastic body is given in the matrix form as 

+ [C]{C/} + [ii:i{C/} = {B} (3.1) 

where [M],[C] and pK] are the mass,damping and stiffness matrices respectively and {R} 
is the external load vector; {U},{?7}and {U} are the displacement, velocity and acceler- 
ation vectors respectively. In the above equation 
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[M] = T,[M\’ 

\K\ = nKr 

{«} = E{B}' 

where the elemental matrices are defined as 

[M]’ = J p [Af]’’ (JV] dV 

m = jT [Bf [D\ [B] dV 

{«}' = / I^^Kr} ds 

In these equations [N] is shape function matrix, [JB] the derivative of shape functions, 
[D] the elastic constitutive relation matrices, and {T} the traction vector. 

In the analysis the damping forces were neglected and therefore the governing equation 
becomes 


[M]{U} + [K]{U} = {i?} . (3.2) 

The above equation can be solved either by time integration or by mode superposition, of 
which the the former is preferred for simulating wave propagation problems of the present 
kind. In this integration scheme, there are many different methods. In the present work 
Newmark integration method for time variable is used. Since the specimen width is large 
in comparison to its thickness ( Fig 2.3 ) plane strain relations are used. Furthermore the 
front face of the specimen remains fixed to the rigid block and its effect on the J-integral 
is negligible (Verma (1995)). The analysis is therefore carried out only on the cantilever 
plate to avoid unnecessary computation. 
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3.2.1 Constitutive relation for Composites 

Verma (1995) developed the program for isotropic materials. It has been modified by Ra- 
makrishna(1997) and Babu(1998) to apply it to the composite materials. In the isotropic 
materials there are only two independent elastic constants, where as in case of the com- 
posite materials there are four independent elastic constants. The elastic constants are 
longitudinal modulus El, transverse modulus Et, rigidity modulus Glt, the major Pois- 
son ratio ult aud the minor Poisson ratio vtl', ojfiy four of them are independent. The 
principal material directions of the composite materials are as shown in Fig 3.1. Tests 
are conducted to find El, Et, Glt and i'lt- El and Et can be directly obtained from 
the tension tests on 0° and 90^^ ply laminates respectively. And one more tension test is 
conducted on a 45'^ ply laminate to get the modulus in that direction. Glt is obtained 
from the expression 

1 1,1 1 1 2,ult 

1^0 ^ A^Wl ^ ^ ^ 

The details of the specimen used and the related results are discussed in detail in Appndix 
A. The obtained values from the experiments are as follows 


El = 40.53 GPa, 

Et = 6.98 GPa, 

Glt ~ 2.84 GPa and 
Ult — 0.277 . 

The constitutive relation for the composite materials is given by 


= [Q] {^} 


where 


38 



[Q] 


Qn Qu 0 
Qi 2 Q 22 0 

0 0 Qee 

Qij axe the stiffness coefficients of the composite. They can be expressed as 


Q 


El 


11 


1 — Ult 


Q 22 — 


Ej 


1 — Uix ptl 


(3.4) 


Q\2 — 


utl El 


vlt Et 


1 — Ult ^TL 1 “ ^LT ^TL 


Qee = ^LT 

3.2.2 Path independent integral 

Under appropriate assumptions of material homogeneity, the strength of the crack tip field 
is governed by an integral evaluated over a path that is far removed from the crack tip. 
Since the stress and displacement data axe evaluated far away from the crack tip, these 
axe relatively insensitive to the finer details of modelling of crack tip region. Therefore 
a path independent integral,(J integral) is used in the analysis. In this analysis the J 
integral given by Kishimoto, Aoki and Sakata (1980) is used. It is defined as 

J = Urn, ^0 / [W - TiUij] ds + [ p Ui dV 

7rc-+r+rc+ 

where W is the strain energy density, Ui is displacement, Uj is the direction cosine of the 
unit outward normal in x\ direction, Tj is the traction, Tc , T, are defined paths 
(Fig 3.2) and V, are voltunes. 
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The value of J is independent of the choice of the path only under ideal conditions 
such as steady state crack growth. In many cases though these conditions are not met, 
the path independence of the J integral can still be established within an allowable error. 
The path was held stationary as the crack tip is extended in a self similar manner. 

3.2.3 Parameters in the Program 

The choice of certain parameters will infleuence the outcome of the program. Such pa- 
rameters in this problem are time step, mesh size, etc. The choice of time step (At) for 
time integration is important for an accurate solution. An optimum choice of time step 
is At = d/c, where d is the smallest mesh size and c is the fastest wave velocity. However 
it has been observed that the the value of J integral does not vary much with time before 
opening the crack. Therefore in the present study, a coarse time step of 0.2^s is selected 
prior to the growth of the crack and fine time step of 0.02 /xs is used after the crack 
growth in initiated. 

Mesh size 

As discussed by Seron et al. (1990) higher order elements are not required in the analysis 
of dynamic problems. Thus four noded isoparametric elements are used in this work. 
To take into the consideration of bending as well, fine mesh is used. At the same time 
to avoid non uniformity of mass distribution which leads to reflection of waves, umform 
mesh is used throughout the specimen. Figure 3.3 shows mesh generated in the present 
work. 


3.3 Crack opening scheme 


As the crack propagation velocity is lower than the wave velocity, the crack tip occupies 
positions in between the nodes at various time steps of the numerical integration. Thus 


if a simple node shifting procedure is used, then the crack tip either stay at one node or 
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jump from one node to the next node in time At in the computational simulation. When 
the crack tip movement is irregular spurious fluctuations are superposed on the value of 
J determined through simulation. 

Thus to overcome this difficulty Verma (1995) adopted a technique where the nodes are 
released gradually. The holding force at any node is linearly reduced to zero as the crack 
reaches the end of the next element. If the actual crack tip is located at C in between 
the element nodes B and D as shown in Fig 3.4. The holding back force F at node B is 
gradually reduced to zero as the crack tip reaches the node E. When the crack tip is in 
between nodes B and D, the holding force is given by, 


— = [1 — —] 
Fbc ^ 2 


(3.5) 


where Fbc is the force at node B when the node was not opened and the crack tip was 
located at point B, hi is the crack extension beyond node B and d is the element length 
as shown in Fig. 3.4. When crack propagates beyond the node D to a point Di 


Ib_ ^ r^_ d + b2 , 

Fbc ^ 2d ^ 

3.4 Initiation and Propagation toughness 


(3.6) 


Initiation fracture toughness is the fracture resistance of the material tmder dynamic 
loading conditions for a specific crack length and specified boimdary conditions. Initia- 
tion fracture toughness represents the toughness of the material just before the initiation 
of crack. Figure 3.5 shows the variation of J integral with time for the stationary crack. 


At the initiation time, determined experimentally, the crack opening scheme is is invoked 
(Eqs. 3.5,3. 6) in the numerical simulation. Figure 3.6 shows rise of J for stationary 
crack upto the initiation time and hatJ for the propagating crack. The toughness at 
the initiation time is called the initiation toughness {Jini)- As shown in Fig 3.6, once the 
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crack is initiated, J drops suddenly to lower values and after some time it stabilizes. This 
stabilized value of the J integral is treated as the propagation toughness ( Jprop )■ The 
results are discussed in detail in Chapter 4. 

3.5 Summary of the program 

Figure 3.7 shows the flow chart of the program developed by Verma(1995). The input for 
the program is end deflection of the cantilever u{t), determined experimentally. At the 
initiation time, crack length as function of time is provided to the computer simulation 
and crack initiation through node release is invoked. 

3.6 Closure 

The aim of this chapter was to describe the finite element formulation, the path indepen- 
dent integral, the crack opening scheme, the initiation and propagation toughness, and 
the summary of the program. 
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Figure 3.1: Principal directions of composite material 
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Figure 3.3: Mesh generated for the simulation 
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Figure 3.4: Crack opening scheme 
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Figure 3.7: Flow chart of the FE code 
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Chapter 4 

Results and Discussion 


4.1 Introduction 

Experiments are conducted to determine the deflection of the cemtilever end, crack ve- 
locity and initiation time. These are the inputs to the FE code developed by Verma 
(1995) along with the material properties. The output of the program is the variation of 
J-integral with time in different conditions such as the stationary crack and the propagat- 
ing crack. This chapter briefly describes the results of initiation toughness, propagation 
toughness and their comparison with the quasistatic interlaminar toughness and dynamic 
toughness determined by other investigators. 

The precrack length , position of the strain gauges as measured with the travelling mi- 
croscope are tabulated in Table 4.1. 

4.2 Interlaminar Initiation Toughness 

Interlaminar initiation toughness of the material is the toughness of the material when 
the crack is about to initiate. It is the fracture resistance of the material under dynamic 
loading conditions for a specific crack length and boundary conditions. This section de- 
scribes the results of the initiation toughness. Experiment 1 is discussed in detail and 
results of other four experiments are briefly presented highlighting the points when ever 
necessary. 
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Table 4.1: Specimen thickness, width, precrack length and strain gauge locations for 
Expt. 1 to 5. 


Expt. No. 

Specimeu 

Thickness 

(mm) 

Specimen 

Width 

(mm) 

Precrack 

Length 

(mm) 

Strain gauge locations 

SGI 

(mm) 

SG2 

(mm) 

1 

3.60 

25.5 

41.4 

43.8 

46.4 

2 

3.20 

25.1 

39.3 

41.6 

43.9 

3 

4.10 

25.9 

40.1 

41.9 

44.6 

4 

4.10 

25.8 

40.3 

42.3 

45.1 

5 

4.04 

25 

39.6 

41.8 

44.5 


central 

Experiment 1 iiftiKlk Mk 1 £ 6 £ J g 

When specimen of 3.6 mm thickness, and precrack length of 41.4 mm is impacted by a 
load bar part of the energy is transmitted to the specimen causing the deflection of the 
cantilever and the remaining energy is reflected back as a tensile pulse. Strain gauges 
mounted on the load bar monitors the stress pulses and the strain gauges bonded ahead 
of the crack tip are used to monitor the crack propagation. The response of the load 
bar strain gauges and and the velocity strain gauges as recorded in the oscilloscope are 
as shown in Fig 4.1, Channel 1 gives the record of the incident and reflected pulses in 
the load bar. Channels 2 and 3 give the response of the first and second velocity strain 
gauges (SGI and SG2) respectively. The velocity of the load bar end is calculated using 
the expression (Sec 2.3,5) 
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V2{t) = 


CTj — (Ti 


P c 

where <ti is the compressive incident pulse and 0-3 is the tensile reflected pulse in the load 
bar; p, the density’ of the material and c, the longitudinal wave velocity. The variation 
of velocity with time is shown in Fig 4.2(a). The end displacement of the load bar is 
obtained from the following expression, which is obtained by integrating the load bar 
velocity with respect to time (Sec 2.3.5). 

Uo{t) =[ Vo dt = [ -) 

Jo Jo p c 

Since the specimen is screwed to the load bar the end displacement of the load bar be- 
comes the deflection of the cantilever end. The variation of deflection of the cantilever 
end with time is as shown in Fig 4.2(b). The loading of the cantilever starts at the time 
corresponding to mid point of the head of incident pulse and head of reflected pulse. This 
is denoted as "REFERENCE" in the Fig 4.1. 

The response of the velocity strain gauges leads to the determination of crack velocity. 
When the crack tip passes near a strain gauge the singular strain field gives a strain peak 
that is recorded in the oscilloscope. Figure 4.3 shows the blown up view of the response 
of the velocity strain gauges. The crack velocity is determined from the time difference 
of the peaks of the strain gauges and the distance between the two strain gauges (Sec 
2.3.6). The strain gauges are 2.6 ttitti apart and the difference in the peak response is 
5.5 ps. Thus the velocity is found out to be 472 m/s. The initiation tune is obtained by 
extrapolating the crack velocity up to the precrack length (Sec 2.3.7). Figure 4.4 shows 
the extrapolation to determine initiation time for Expt.l. The first strain gauges in this 
experiment is located at a distance of 2.4 mm from the precrack. Thus the initiation 
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time obtained is 50.02 fjs. 


These experimental data (the deflection of the cantilever end. crack velocity and initi- 
ation time) axe used as input to the FE code as discussed in chapter 3. The variation 
of J integral in the stationary phase is as shown in the Fig 4.5. The value of J builds 
up and then drops down. The rise of the interlaminair toughness till the initiation is also 
obtained experimentally by John Lambros and Rosakis A.J (1997a). They investigated 
the variation of dynamic fracture toughness up to the initiation time on unidirectional 
graphite/epoxy composite plates with edge notch under impact loading (Fig 4.6). In the 
present study the experiment is controlled in such a way that the initiation occurs before 
the J reaches the peak value. in this experiment is found out to be 149 J/m* (Fig 4.5). 

Experiments 2-5 

Similar to the Expt 1 four more experiments were carried out to determine the initiation 
and propagation toughness. The strain gauge peaks, crack velocity initiation time and 
Jini for various experiments are tabulated in Table 4.2. 

Experiment 2 is carried out on a relatively t hinn er specimen (3.2m77i) and with a precrack 
length of 39.3 mm. Figures 4.7-4.9 show the oscilloscope traces, variation of velocity of 
load bar with time, deflection of cantilever end with time, and the details of extrapolation 
for this experiment. The crack velocity is 657 m/s and the initiation time is 51.5 ^ls, 
Figure 4.10 show the variation of J integral with time. The initiation toughness {Jini) is 
154 J/m 2. 

Experiment 3 is carried out at high crack velocity. The oscilloscope traces and the 
processed results are shown in Fig 4.11- 4.13. The crack velocity is 900 m/s and the 


Table 4.2: Crack velocin-. initiation time and for different experiments. 


Expt. No. 

Precrack 

Length 

(mm) 

SGI from 
Cracktip 
(nmi) 

Peak Re.sponse 

Crack 

Velocity 

{m/s) 

Initiation 

Time 

il^s) 

Jini 

(J/m^) 

SGI 

ins) 

SG2 

(ns) 

1 

41.4 

2.4 

54 

59.5 

472 

50.02 

149 

2 

39.3 

2.3 

55 

58.5 

657 

51.50 

154 

3 

40.1 

1.8 

52 

55.0 

900 

50.00 

357 

4 

40.3 

2.0 

51 

54.0 

933 

48.86 

322 

5 

39.6 

2.2 

55 

58.0 

900 

52.42 

141 


initiation time is 50 fxs. Figure 4.14 shows the variation of J with time. The initiation 
toughness, being 357 J/nt^ is highest among all experiments of this study. 

Experiment 4 is very much similar to Expt.3. It is conducted at a velocity of 933 m/s. 
The initiation time is the lowest among all the experiments. It’s value is 48.86 fis. Figures 
4.15-4.17 shows the traces and processed data. The variation of J integral with time is 
shown in Fig 4.18. The initiation toughness is 322 J/m^. 

Experiment 5 is also carried out at high crack velocity. But unlike in the case of previous 
two experiments where the velocity is high and the initiation toughness is also high, the 
initiation toughness is only 141 J/m?. Figures 4.19-4.21 shows the oscilloscope traces 
and processed data. The variation of J integral is shown in Fig 4.22. 
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4.3 Interlaminar Propagation toughness 

For determining the interlaminar propagation toughness (Jprop) the input to the finite ele- 
ment code consists of the crack propagation history apeirt from the variation of cantilever 
end deflection with time. While executing the FE code the crack propagation algorithm 
is called after the initiation time. Similar to the procedure of the previous section, Expt.l 
is discussed in detail and other results are presented briefly. 

Experiment 1 

As discussed by Verma(1995), once the FE code is executed upto the initiation time, the 
crack propagation module is called. The location of the crack tip in each iteration is 
known through the experiment; an appropriate factor is applied to the holding back force 
for modelling crack propagation (Eqs. 3. 5, 3. 6). A tune step of 0.02/iS is chosen in the 
propagation phase. Using the boundary condition and the crack propagation data the 
dynamic FE analysis deter min es the stress/strain field in the specimen in successive time 
steps. Then the variation of J integral with time is obtained as shown in Fig 4.23. In the 
be gin in g the crack remain stationarj' upto the initiation time and J integral increases. At 
the initiation time the crack tip starts growing under the known dynamic displacement 
boundary conditions. Figure 4.24 shows the blown up view of the variation of J integral 
in the crack propagation stage. The J integral stabilises after sometime and the stabiUsed 
value is treated as the propagation toughness (o^op). The propagation toughness in this 
case is 70 J/m^. 

The oscillatory behaviour of J integral is due to several reasons (Verma (1995)). First 
and the foremost, the free surfaces of the specimen are very close to the crack tip and, 
in fact, the propagation time for some stress waves to emanate from the crack tip and to 
return to the crack tip are as small as 2 - 3^s. The superposition of these waves and their 
effect on J integral will provide ripples. Second, in the finite element modelling the crack 
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Table 4.3: Values of crack velocity and Jp^qp for different experiments 


Expt.No 

Crack velocity 

77i /s 

Propagation Toughness(Jj„.qp) 
J/m^ 

1 

472 

70 

2 

657 

42 

3 

900 

80 

4 

933 

95 

5 

900 

30 


tip is moved from one point to another point in discrete steps by linearly decreasing 
the holding back force. Thus the FE code may not be simulating the crack propaga- 
tion behaviour very accurately. The nodal release mechanism can further be modified to 
smoothly release the crack to rule out its effect on the oscillatory nature of J^rop- 

Experiments 2-5 

Variation of J integral with time for the propagating crack for Expt.2 is shown in Fig 
4.25. Figure 4.26 shows the blown up view of the J variation in the propagation phase. 
Jprop for this experiment is A2Jlm"^. Table 4.3 gives the values of the Jprop for different 
experiments. 

In experiment 3, where the crack velocity and the initiation toughness are high, the prop- 
agation toughness is also high. Figure 4.27 shows the variation of J of the dynamic crack. 
The blown up view of the variation of J in propagation phase is shown in Fig 4.28. 
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The response of Expt.4 is similar to Expt.3. The value of Jprop here is 95 J/rn^, the 
highest among all the values. The variation of J and its blown up views are shown in 
Figures 4.29, 4.30 respectively. 

Figure 4.31 shows the variation of J for Expt.5. The blownup view of the same is shown 
in Fig 4.32. Even though the velocity is high in this case, the value of J^-cp is lowest 
among all the values ( 30J/m‘). 

4.4 Comparison of Dynamic interlaminar toughness 
with quasistatic interlaminar toughness 

To compare the dynamic interlaminar toughness with quasi static interlaminar toughness 
experiments are conducted to determine the quasistatic interlaminar toughness {Gjc) 
through a standard test (Guedra.D et, al (1987)). The tests involves pulling a DCB 
specimen in displacement control mode. Load is applied till the crack extends by a small 
distance ( 5-10 mm). Then the machine is stopped till the crack becomes stationary 
and then the specimen is unloaded. The specimen, now with a longer crack length, 
is loaded again. The compliance of the specimen is determined through the loading 
curve. The specimen is subjected to several loading-unloading cycles so as to determine 
the compliance and critical load for several crack lengths. With proper data reduction 
scheme (Guedra.D et. al (1987)), the critical energy release rate is evaluated. Value.s of 
Gjc obtained from the experiments are tabulated in Table 4.4. 

It is clear from the Tables 4.2, 4.3, 4.4 that the quasistatic interlaminar toughness is 
significantly higher than the dynamic interlaminar toughness. Thus it is important to 
characterize the interlaminar fracture toughness in dynamic conditions. 
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Table 4.4: Experimentally obtained values of Gjc 


Expt.No 

Quasistatic interlaminar 
Toughness (G/c) 

J/m^ 

Average value 1 

1 

1304 


2 ' 

1020 

1059 

3 

790 


4 

1124 



4.5 Comparison of Results 

In this section the values of initiation toughness and propagation toughness are compared 
with those of Verma (1995), Ramakrishna(1997) and Babu(1998). 

The variation of Jini and J^ap for the present study is shown in Fig. 4.33. The quasistatic 
toughness is much higher than the dynamic interlaminar toughness and also imtiation 
toughness is substantially higher than the propagation toughness. 

Figure 4.34 shows the variation of imtiation and propagation toughness as obtained by 
Verma (1995) on a specimen made of two steel strips bonded together with epoxy. It 
is worth noting that decreases with increasing velocity of interlaminar crack. Also 

It is interesting to observe that unlike the present study, the quasistatic toughness is 
significantly lower than the initiation and propagation toughness. It may be due to the 
weak adhesion between epoxy and steel surface. In the present study the resin which 
binds the plies together is same as the resin of the bulk material, thus heterogeneity is 
of much smaller degree. Thus micromechanisms of failure are different in nature and 
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crack growth in steel/epoxy specimen has a threshold level that is much higher than the 
quasistatic toughne-ss. However, even in steel specimen Jprc^ is substantially lower than 
the quasistatic toughness. 

The \-ariation of J with crack velocity as obtained by Ramakrishna (1997) is shown m 
4-35- lit his work GFRP laminates used were not as superior as the one used in the 
present study. The quasistatic toughness ( 415 J/m") is higher than the initiation and 
propagation toughness. However each of these value is lower than those of the present 

Study. 


The results obtained by Babu(1998) are as shown in Fig 4.36. The initiation toughness 
was 245-600 J/m^ and the propagation toughness was 45-83 J/m-. These values are very 
much comparable to those of the present study, where the initiation toughness is 140-360 
J/m' and the propagation toughness is 30-95 J/m'. Figure 4.37 shows the comparison 
of J,,. for the present work with that of Babu’s work and Fig 4.38 shows the comparison 
of for same. The low value of the initiation toughness in the present study may be 
due to the high crack velocity. 

4.6 Closure 

The aim of this chapter was to find the initiation toughness and propagation toughness 
of interlaminar cracks in GFRP laminates. The experimental measurements such as the 
end deflection of the cantilever, crack velocity and the initiation time are used as the 
input to the FE code to simulate the dynamic fracture behaviour of the GFRP laminate. 
The initiation toughness was found to be 140-350 J/m' and the propagation toughness 
is 30-95 J/m'. To compare the results with the quasi static interlaminar toughness tests 
are conducted to find the same. It was found that the dynamic interlaminar toughness 
of the laminate is much smaller than the quasi static interlaminar toughness. 
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figure 4.1; Oscilloscope Traces of Expt. 1 
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ca) 



(b) 

Figure 4.2: (a)Variation of velocity of load-bar-end with time for Expt. 

(b) deflection of cantilever end with time for Expt. 1 


1, and 
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Figure 4.3; Blown up view of responses of Velocity strain gauges( Expt.l) 
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Figttre 4.7: Oscilloscope Traces of Expt. 2 
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Figure 4.11: Oscilloscope Traces of Expt. 3 
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Figure 4.15: Oscilloscope Traces of Expt. 4 
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Figure 4.19: Oscilloscope Traces of Expt. 5 
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Figure' 4.3.H: Variation of J witii crack velocity for the present work 
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Figure 4,34: Variiitiou of J with crack v<'locity, as obtained by Verma (1995) 
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Figure* 4.3r>: Variation of ./ with crack v<‘lo<’it.y as obtairuni by Ramakrishna (19.)7) 
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Figure' 4.38: Comparisiou of Propagation toughness with that of Babu(1998) 
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Chapter 5 

Conclusions and scope for further 
work 

5.1 Conclusions 

A (Icwclojx'd by Verina(1995) has been applied to determine the initiation and 

propagation toughn<>ss of fast moving cracks in GFRP laminates. This is a combined 
experim(mtal and numerical technique. 

The experimental technique consists of a specimen; the front face of which is bonded to 
a rigid block and the cantlever end on the rear is screwed to a load bar. A striker bar 
impacts th(’ load bar, wliich in turn impacts the specimen and causes the crack propaga- 
tion. Two strain gauges bonded on the load bar monitors the stress pulses, which lead to 
the d('tennination of th<' end displacement of the load bar. Since the specimen is screwed 
to the load l)ar, the end chsplacernent of the load bar becomes the deflection of the can- 
tik'ver end. Two strain gauges bonded ahead of the crack tip on the specimen monitors 
the crack propagation. Velocity of the crack is determined from the strain peaks in the 
respons(^ of these strain gauges and the distance between them. It is then extrapolated 
to the precrack length to determine the crack initiation time. 
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Thus the experimentation provides the deflection of the cantilever end, crack velocity and 
the initiation time. These along with the material properties are used as input to the FE 
code to simulate' the fast moving crack. Gradual nodal releeise method is used to model 
the crack propagation. The initiation toughness ( Ji„i) was found to be 140-350 J/m^ and 
the propagation toughness {Jprop) to be 30-95 J /rn?. These results are then compared 
with quasistatic interlaminar toughness ( G/c). It has been observed that the dynamic 
interlaminar toughness of the crack is much smaller than the quasistatic interlaminar 
toughness. 

5.2 Scope for further work 

• Insti'ad of strain gauges thin conducting strips can be used ahead of the crack tip 
to monitor tlu' crack velocity 

• Th(' precrac'k length and the thickness of the cantilever can be varied to find their 
effect on the values of Ji„i and Jprop- 

• The work may he extended to angel ply laminates. 
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Appendix A 


For an orthotropic composite laminate the principal material directions are the longitu- 
dinal direction ( along the fibres ) and transverse direction ( across the fibres ). Unlike 
the case of isotropic materials where there are only two independent material constants, 
in this case there are four independent constants. The constants are the longitudinal 
madulas ( Ei), transverse modulus {Ej), shear modulus [Git) and the major poissons 
ratio {vlt). Ei and Ep are determined by conducting tensile tests on 0° and 90“ laminates 
respectively. u,;r can be determined from the tensile test on 0“ laminate by measuring 
the strain in longitudinal and transverse directions. Figure A.l shows a sample stress- 
strain curve for a tensile test on 0“ laminate. Figure A.2 shows the same on 90“ laminate. 
Tensih' test conducted on 45“ laminate yeilds E,^. Git can be found firom E45 using the 
formula, 

^ i/JL 4. _1_ _|_ ^ 

4 El Ep Git El ^ 

Figure A. 5 sliows a sample stress-strain diagram for tensile test on 45“ laminate. 


The tests ar<' conducted on 16 ply laminates. The values olEi, Ep , Gip and vip obtained 
from the experiments are tabulated in Tables A.l, A.2, A. 3, and A. 4 respectively. 


Tal)l(' A.l Exja'rimental vahies of Ep 


S.No 

Longitudinal Modulus [Ep) 
{GPa) 

1 

46.30 

2 

35.40 

3 

43.32 

4 

37.90 

5 

42.24 

6 

40.07 

7 

38.50 
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Table A.2 Experimental values of Et 


S.No 

IVansverse Modulus {Et) 
{GPa) 

1 

6.35 

2 

8.17 

3 

7.10 

4 

5.80 

5 

7.60 

6 

7.16 

7 

7.40 

8 

7.20 

9 

6.06 


Tabk' A. 3 Tabulatc'd values of G^t 


S.No 

E 45 

{GPa) 

Modulus of Rigidity Git 
{GP a) 

1 

8.20 

2.997 

2 

8.02 

2.902 

3 

7.82 

2.798 

4 

7.50 

2.673 


Tal)k' A. 4 Tal)ulat(‘(l values of iz/r 


S.No 

Major Poissons ratio {i'it) 

1 

0.281 

2 

0.265 

3 

0.287 
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Figure A. I : Stress-strain diagram for tensile test on 0“ laminate 
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Figure A. 2: Stress-strain diagram for tensile test on 90'" laminate 
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Figure A.3 : Stress-strain diagram for tensile test on 45° 


laminate 
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